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DARK FLUID: TOWARDS A UNIFICATION OF EMPIRICAL THEORIES OF GALAXY ROTATION, 
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ABSTRACT 

Empirical theories of Dark Matter like MOND gravity and of Dark Energy like f(R) gravity were 
motivated by astronomical data. But could these theories be branches rooted from a more general 
hence natural framework? Here we propose the natural Lagrangian of such a framework based on 
simple dimensional analysis and co-variant symmetry requirements, and explore various outcomes 
in a top-down fashion. Our framework preserves the co-variant formulation of GR, but allows the 
expanding physical metric be bent by a single new species of Dark Fluid flowing in space-time. Its 
non-uniform stress tensor and current vector are simply functions of a vector field of variable norm, 
resembling the 4- vector electromagnetic potential description for the photon fluid, but is dark (e.g., by 
very early decoupling from the baryon-radiation fluid). The Dark Fluid framework naturally branches 
into a continuous spectrum of theories with Dark Energy and Dark Matter effects, including the f{R) 
gravity, TeVeS-like theories, Einstein-Aether and i^A theories as limiting cases. When the vector field 
degenerates into a pure Higgs-like scalar field, we obtain the physics for inflaton and quintessence. 
In this broad setting we emphasize the non-constant dynamical field behind the cosmological con- 
stant effect, and highlight plausible corrections beyond the classical MOND predictions. Choices of 
parameters can be made to pass BBN, PPN, and causality constraints. The Dark Fluid is inspired to 
unify/simplify the astronomically successful ingredients of previous constructions: the desired effects 
of inflaton plus quintessence plus Cold DM particle fields or MOND-like scalar field(s) are shown 
largely achievable by one vector field only. 

Subject headings: Gravitation; Cosmology: theories; Dark Matter; Galaxies: kinematics and dynamics 



; 1. INTRODUCTION 

OO ' Gravity, the earliest and the weakest of the known 
00 . forces, has never been very settled. The beauty of co- 
lO variant symmetry motivated Einstein to supersede the 
' Newton's paradigm with General Relativity, but (inade- 
, quate) empirical evidences motivated Einstein to intro- 
■ duce first and then withdraw the cosmological constant, 
a concept defying quantum physics understanding even 
. in modern day. While making generally tiny or a factor 
;V ' of two corrections, the equivalence principles insist on 
.5^ I certain symmetries in space-time, e.g., co-variant sym- 
, metry and no frame to measure locally any absolute di- 
^ rection of gravitational acceleration for a free-falling ob- 
server. However, symmetry can be spontaneously broken 
if there are dynamical interactions or couplings of fields; a 
well-known mechanism in several branches of physics, es- 
pecially the Higgs-mechanism in particle physics to give 
a mass to a particle. Many attempts have been made 
to break the strong equivalence principle by adding new 
fields (degrees of freedom) in the gravity sector, which 
essentially means the gravitational "constant" G may 
be is a new dynamical degree of freedom governed by 
other fields coupled to the metric. The best known is 
the Brans-Dicke theory (1961). The lesser known is that 
a vector field of non-zero absolute value in vacuum can 
also be coupled to gravity, to give absolute directions 
(Will 1993). It has long been suggested that Lorentz 
symmetry can be broken locally in the quantum gravity 
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and string theory context (Kostelecky & Samuel 1989) to 
yield a vector field of a non-zero expectation value {e.g., 
pointing towards the direction of time) in the vacuum. 
The most successful attempt so far is the Einstein- Aether 
theory of Jacobson et al. (2001). A common theme 
of these theories is that they are not invented for cer- 
tain observational anomaly. Rather in the same vein as 
how symmetry motivated General Relativity, these the- 
ories meet the astronomical data only a posteri, e.g., Li, 
Mota, Barrow (2007) showed a vector field in the gravity 
sector could NOT be excluded by the accurate Cosmic 
Microwave Background (CMB) data. 

Nevertheless, the above order is not the only way to 
discover theories. The puzzling black body radiation 
spectrum and Balmer's curious empirical formula for hy- 
drogen lines are among the odd pieces of classical physics, 
which lead to full formulation of quantum mechanics. 
This bottom-up approach is often gradual, the arrival 
of the final theory taking several generations of formu- 
lations {e.g., from Planck's model for blackbody radia- 
tion and Bohr's model for hydrogen atom to Heisenberg's 
matrices-based formulation in general) with different lev- 
els of mathematical rigour and sophistication. 

Historically, Milgrom's Modified Newtonian Dynam- 
ics (MOND) was invented without any packaging by co- 
variant theories of gravity, just as the Dark Matter em- 
pirical concept was invented by Zwicky without pack- 
aging first with SuperSymmetry-like particle field theo- 
ries. MOND is a formula motivated to reveal the curi- 
ous uniform rules (or facts) underlying rotation curves 
of many spiral galaxies, as Balmer's formula and its gen- 
eralisations suggesting strongly a fundamental rule for 
all atomic lines. Since the rule is empirical and bare 
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(without CO- variance) , it waits to be enshrouded by a 
theory preserving basic symmetries to predict any log- 
ical corrections to situations where the empirical rule 
must fail slightly, e.g., by factor of two in some gravi- 
tationally lenses made by elliptical galaxies and clusters 
of galaxies. The Tensor- Vector-Scalar (TeVeS) frame- 
work of Bekenstcin (2004), unifying earlier constructions 
by Sanders and others, makes the first step to the in- 
tegration of MOND formula with fundamental physics. 
A time-like vector field is shown to be the necessary in- 
gredient of a MOND gravity. Yet the original aim of 
TeVeS was limited, e.g., not addressing the cosmological 
constant problem, or the inflation. Orthogonally many 
literatures considered theories of modified gravity such as 
the f{R) gravity (Chiba 2003) and scalar inflation the- 
ory as ad hoc fixes of the cosmological constant problem 
and the horizon problem respectively, without aiming to 
address outstanding questions on galaxy rotation curves. 
Recently Zhao (2007), built on the work of Zlosnik et al. 
(2007), showed that these outstanding problems of DM 
and DE can find at least one common solution simulta- 
neously in the framework of a massive vector field. In 
these theories, there is "One Field which rules them all 
and in the darkness bind them." 

The most famous examples of a vector field is the mass- 
less spin-1 photon and the massive Z-boson in the electro- 
weak theory. The standard way to give masses to parti- 
cles in particle physics is the Higgs mechanism where a 
scalar field, coupled to the vector field, acquires a non- 
zero value in vacuum. Turning to the gravity sector, 
however, it is unclear how the fundamental prediction of 
spontaneous symmetry breaking from quantum gravity 
might be related to the mundane effects of galactic dark 
matter and cosmological constant. We lack a framework. 

We shall show that MOND, /(i?) gravity, Einstein- 
Mther theory and inflation can be integrated into a com- 
mon framework with a unit vector and a dynamical scalar 
field. MOND would become a specific choice of the po- 
tential of the scalar field. Having such a framework al- 
lows one to explore the consequences of modified gravity 
systematically. It can be meaningless to even differenti- 
ate dark energy and modified gravity. Modified gravity 
contains extra fields, which can be treated as dark energy 
field. 

The goal of this paper is to show the existence of a 
very general Lagrangian for which the MOND formulae 

are the natural consequences in spiral galaxies in equilib- 
rium, rather than the golden rule for (non)-equilibrium 
systems of all scales. We demonstrate this with the modi- 
fied Poisson equation and for the equation for the Hubble 
expansion. 

The outline is as follows. We propose our general La- 
grangian in §2, and illustrate how it reduces to vari- 
ous special cases, TeVeS, BSTV, Einstein-Aether, ,f(K), 
f{R), infiation. We choose a subset of models with 
MOND and Dark Energy effects in §3. We examplify 
the properties of our dark fluid in the case of Hubble ex- 
pansion and inflation (§4), and for static galaxies (§5). 
We discuss corrections to MOND in §6, and summarize 
the properties of the Dark Fluid in §7. Appendix gives 
an estimate of the damping frequency of the Dark Fluid. 

2. THE PROPOSED LAGRANGIAN FOR THE 
DARK FLUID 



Denote a vector field by Z", which has generally a vari- 
able or dynamic norm 

= g,,Z''Z\ (1) 

which is essentially an auxiliary scalar field character- 
izing the norm of the vector field Z", hence is not an 

independent dynamical freedom. A generic coupling of 
the vector field Z° with the space-time is through the 
contractions among the Z°'Z'' tensor, the gab metric ten- 
sor, the Ricci tensor Rat- Hence the most generic theory 
of the vector should contain the terms 

L = [g"b + g-bfi^-" - hZ'^Z'] R^t + hi, (2) 

/34 = /3<^2 + /4VWa<^ + --- (3) 

where fi could be constants or functions of ip. The term 
f2Z'^ZbRab = /2^" ("VaVc^" + VcVa^") could always 
be recast through a full divergence to —^ K^I^S/aZ''VcZ'' 
where K^ji^ are constant tensors or tensor functions. The 
eqiiivalent expression for the Lagrangian is 

L = K^lyaZ''ycZ'' + OoR + 61(^2 ^ e^vVVa^^ + •(4) 

where R = g'^^Rab is the Ricci scalar, and K^jj^ are con- 
stant tensors for simplicity but can also be lengthy func- 
tions of if in general; Oi are constants in the simplest 
case but functions of ip in general. 

The above Lagrangian is generic enough, and many 
dark energy models can be derived from it. For exam- 
ple, the terms (1 + (p'^)R + uj(p'^^^^'^"'> can lead to an 
R + const/i?" (Li & Barrow 2007) gravity (as could be 
checked by solving (p from the equation of motion of ip 
and then substituting back into the original Lagrangian); 
and the terms 4'R — 4(1 — (j)~^)ujV^(pVa'P leads to the 
Brans-Dickc theory of gravity, where the auxilliary field 
4> = l + (p'^. The usual inflation-like or quintessence theo- 
ries can be recovered from the terms R+V{(p)+'^°'(pVa'p. 

The essential dynamics of a cosmological vector field is 
described by the first term KfjiV aZ^V cZ'^ (Ferreira et al. 
2007, Halle, Zhao & Li 2008). If insisting that the field 
has a unit norm guaranteed by a Lagrange multiplier, 
one would recover Einstein-iEther theory (Jacobson et 
al. 2001) and its generalizations (Zlosnik et al. 2007, Li 
et al. 2007). 

3. MOND-INSPIRED SUBSET OF DARK FLUID 
THEORIES 

To see the relation to MOND, consider the subset 
where ©o = 1, ©i ~ Aq, and ©2 ~ A''"^. Let's de- 
compose the foiir dynamical freedoms in the vector field 
into a unit norm part of 3 degrees of dynamical freedom, 

M^ = —, \ = ^^=gabZ''Z\ (5) 
f 

and rewrite the scalar field ip in term of the new scalar 
field A. The Lagrangian is then casted to a form contain- 
ing at least the following 

C = L^ + R + L^ + Lm (6) 

= -^^V||^V||^ - ^ iya^) (VV) + 2AoF(A) + • (7) 

= C4V||^cV,l^'= + C2 (Va^")^ + (^"^a " 1)^* + ' ' ' , (8) 

where N, F, c^, C2 and C4 axe various coupling constants 
or functions of A = ip'^, and Vy = jE"Va, V" = g'^'Vc, 
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S^a — gac-^'^, R is the Ricci scalar, and L* is the La- 
grange multiplier (a kind of potential) . The Lagrangians 
Lm, L^, LjE are for the matter, the scalar field (p and 
the unit vector field jE" respectively, where we omitted 
two possible terms ciKi + C3K3 = ci(VajEh)(V''jE'') + 
C3(Va^fc)(V''jE°) in LjE- Aq is the only dimensional scale 
in the dark fluid, it is a scale of energy density. 

In this Lagrangian we have three dynamical fields: the 
scalar field A, the jEther field jE° and the metric field 
gab pjyg non-dynamical L* . Now varying the action 
^ — ~ J \/—9d'^x 1 itra '^ith respect to them will lead 



167rG 

to the scalar field equation of motion (EOM), jEther 
field EOM and the modified Energy-momentum tensor 
plus the unit vector constraints for M field respectively. 
The general results arc more tedious and are presented 
elsewhere (Halle, Zhao, Li 2008). Here we illustrate the 
physics by considering only the main terms for a specific 
choice of functions. 

3.1. Choices of coupling constants 

The dimensionless function F{X) has the meaning of 
the potential of the scalar field A = <^^, and c^, C2 and C4 
are of order unity, and are also generally functions of A 
and will be shown to be related the sound speeds of the 
Dark Fluid. 



We set the scalar field sound speed 

and we shall treat A'' and as constants. 
We choose the coefficients 



C4(A) 



= A = ^2 



C2(A) _ 6(A) - 1 



1 

"3A' 



(9) 

(10) 
(11) 



i.e. 6=1 — 9?"^. For simplicity we set two other 
terms in of Jacobson's unit vector field to zero, 
i.e., ciKi = and c^Ky, = 0. This might not be 
necessary, but simplifies the analysis of PPN pa- 
rameters in the solar system and the sound speed 
of the vector field. Our choice of Lagrangian with 
Ci = C3 = kills spin-1 mode waves of the vector 
field, and gurantccs that the normal gravitational 
wave in the tensor mode will propagate with the 
normal speed (of light). This choice is perhaps not 
necessary, but is intended to avoid contraversy on 
the causuality issue. Even the spin-0 mode sound 
speed CyE is plausible: a rigourous analysis (Foster 
& Jacobson 2006, valid for any constant A and 6) 
predicts 



, _ C2 (2 



ci) (1-6)(A-1) 1 



C4 (2 -t- 3C2) 



Interestingly in the solar system, where X = (p"^ 
0, the spin-0 mode of the vector field propagates 
almost instantaneously with cje being (SX)^^^"^ 
times bigger than the speed of light, avoiding the 
Cherenkov radiations constraint in the solar sys- 
tem. All PPN parameters are expect to be equal 
to that of GR in the solar system as well; although 
the PPN parameters ai = — 8A, a2 = (3A — 1)A 



arc non-zero, as shown later the scalar field A is ex- 
pected to settle to a very small equilibrium value 
~ (10~^°)^ on earth for our choice of the penalizing 
scalar field potential F{X). 

Our choice for the dark fluid potential (shown in 
Fig.l) is 



Aoi^(A)U=^2 = 



-If 



An 



(13) 



The dark fluid's energy scale is defined by Aq. An 
important property is that in the limit A = — > 1, 

we have 



F' 



d_ 
dX 



F~(l - A)^ if A = (^^ ^ 1, 
ocA-l/^ if A = y.^ ^0; 



(14) 
(15) 

note a prime always means ^). We shall show that 
this property describes a non-uniform (dark en- 
ergy) fluid which gives the MOND-like (dark mat- 
ter) effects in galaxies but not in the solar system. 

4. BACKGROUND COSMOLOGY 

Consider background cosmology in the FRW flat met- 
ric, 

rfs^ = rft2 _ a2(t)(da;2 + dy'^ + dz"^). (16) 

Firstly, the scalar field follows an equation of motion ex- 
actly as a quintessence, 

ip + 3Hip = -{AoF' + 3b'H'^){2N^ip), (17) 

so (p, the norm of the vector field Z'\ tracks the Hubble 
rate H = a/a. The vector field equation of motion gives 
the Lagrange multiplier (or the mass of the vector field) 

1 — 

L* = dt{aH) -\- —^ip^, where a = 26 - 2. This mass is 
varying with time, hence the vector field describes effec- 
tively an unstable slowly decaying particle. 

The modified 00-term of the Einstein equation be- 
comes 

2,H^ = ^ttG[p + pdm+Pde\, (18) 
pDM = {b-^-l)p (19) 



PDE '■■ 



1 



87rG6 



1 



2iV2 



'P^ + AqF 



(20) 



where p is the (background) energy density of baryon- 
radiation fluid. For our choice of 6, 



6-1 = 



p^ 



r 



(21) 



so we get a dark-matter-like effect of the vector field for 
1 < A = </?2 < 2 by amplifying the gravitational constant 
G by a factor 00 > 6~^ > 2. There is no dark matter- like 
effect at very high redshift, e.g., radiation era or BBN, 
where A ^ 00, hence 6^1. So the BBN constraint 
is automatically satisfied because the Hubble expansion 
rate at BBN is equal to that of a radiation only universe. 

A more general choice of dark fluid potential of this property 
is F oc / (i^" — for n = 1 (as above) and n = 2, 3, 4, .... These 

potentials are always simple polynomial functions of ip. 
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Equivalently the Einstein equation can be written as 

(22) 
(23) 



- {2-+H']=8TrG{p+pDM+PDE) 



PDM = {b ^ 

1 

Pde- 



SirGb 



l)p 
1 

2Ar2 



AqF + 4b'Hipip 



?4) 



where the pressure of the baryon-radiation fluid p = in 
the matter-dominated era, hence Pdm = 0. Apply the 
slow-roU approximation in the late universe we find the 
effective pressure of the vector field 



■ Pde ~ Pde 



Ao 



F 



Ao (v 



1 l-' -2 



SnGb SttG ((^ + 1) 



(25) 



This behaves like a dark energy with 

w = Pde /pde ~ -1 (26) 

and with a characteristic scale Aq, which must be set of 
order (8 x 10~^°m/s^)^ to match the observed cosmo- 
logical constant. Note that in writing above equations 
we have implicitly assumed that the effective dark mat- 
ter and effective dark energy components couple to each 
other. This can been seen by checking that neither Pdm 
nor Pde satisfies the conservation law p + 3H{p+p) = 0, 
but their sum does. The coupling strength is determined 
by 6': if 6 is a constant (equals to unity for our model in 
the early universe), then the two components decouple. 

Models with a very large A'' could even drive inflation in 
the early universe when the energy density is dominated 
by the norm of the vector field {i.e., the scalar field (p). 
It is easy to verify the solution 



A, 



3H% 



(27) 



3 dlna 

applies in the slow-roll phase, where 6^1 for very 
large ip. This phase of slow rolling can infiate the uni- 

(2_ 2 \ 
"^gjya^ j (Kanno & Soda 2006). 

For the universe to inflate by a factor exp(60), e.g., tp 
rolles from an initial value (pi ^ 20N to a final value 
~ lOA^. The end of inflation or the start of the radia- 
tion era is at a time ~ (^)"^^^ ~ 10i^F-V2sec~ 

10^ (t^) the time scale can be even shorter for 

other forms of Dark Fluid potential The inflation might 
end when the vector field decays partly into known parti- 
cles via some small coupling perhaps of the type gabZ^\J^ 
between the vector field and the current field of some 
known fields, e.g., coupling with sterile neutrinos, which 
would then mix to neutrinos of all fiavors and couple to 
photons, leptons and hardrons etc. 

5. STATIC GALAXY LIMIT 

To work out perturbations in static galaxies, remember 
that in the Newtonian gauge we have only two scalar 
mode perturbation potentials, $ and \E', which appear in 
the perturbed metric: 

ds^ = (1 + 2^)dt^ - al{l + 2^){dx^ + dy^ + dz^){2S) 



10" 



-9/2 



l)^difi oc (fi^ then R-^ 



wc will let ao = 1- We assume NO Hubble expansion. 

The vector field eqiiation of motion in static systems 
fixes M"- = (1 — 0, 0, 0), so the vector field tracks the 
metric exactly without any freedom in static galaxies. 

The 00-component of the Einstein equation becomes a 
Poisson equation 

^ -{2^),U = ^T^G{p + PDM + bpDE), (29) 



pDM : 



87rG' 



(30) 



where we use notation F^i = diF, and the dummy in- 
dex implies co-variant or contra- variant derivatives with 
respect to x, y, z. Wc use the approximation that the 

■ 2 

DE part SnGbpDE = ^ -|- Aoi^ is a negligible source 
compared to SnGp from the baryons, and that 

- * = $ (31) 

from the spatial cross term of the Einstein equation. The 
above result is essentially a Poisson equation where the 
vector field creates an effective dark matter-like source 
term pdm- Rearrange the terms, the same equation be- 
comes the MOND Poisson equation 

V • [(1 - A)V$] = 47rGp, A = <^^ (32) 

To see that 1 — A can be identified with the MOND 
Pm function, first we define a value of the scalar field 
fM such that 

' ' (33) 



A, 







We find that the scalar field equation of motion is given 
as 

- 4VV = - [AoF' - |V$p] (2ArV), (34) 

where wc neglect all time-dependent terms. This equa- 
tion is similiar to the equation of Yukawa potential with 
a screening length of L, (V^ — L~'^)(p = 0. In the sim- 
plest case, we adopt to kill the Laplacian term 
= y',_^ „ , didi. In the static limit, we find the equa- 
tion for the scalar field becomes 



4Ao(A-V4 - xV^f 



|V$|2 (35) 

for our choice of F{\). The equation can then be solved 
as 

-4 



A: 



2 




(36) 



To see we recover the properties of MOND function 
Pm (see 1 — ip\[ vs x shown in Fig.l), we rewrite the 
solution of the scalar field as 



1 2 _ 

~ = Pm 



where x = 



/a; 



[37) 



l-(f) ,where|V$|» 
This is exactly the physics of MOND if 

Va^ 



ao (38) 

is identified with the MOND acceleration scale oq. In 
the solar system or strong gravity regime, the modifica- 
tion factor 1 — {x/A)~^ ^ 1 to the Newtonian Poisson 
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equation is small and reduces sharply. In weak grav- 
ity, applying spherical approximation around a dwarf 
galaxy of mass mb, we have |V$p/ao = Gmbr~^, and 
the rotation curve V^^^{r) = rV$. The big success of 
MOND in dwarf spiral galaxies is to explain their Tully- 
Fisher relation VX-{r)/{Gmb) = oq ~ IQ-'^^m/s^ if 
Ao ^ (1 X lO^^'^m/s'^)'^ , which is of the order of mag- 
nitude of the observed amplitude of "the cosmological 
constant" effect. In the intermediate regime, our fipj re- 



sembles the "standard" /i 



function of MOND, 



SO it will fit rotation curves of galaxies very well. 

6. TEMPORAL AND SPATIAL CORRECTIONS TO 
MOND: OSCILLATIONS AND DIFFUSIONS 

When considering merging systems like galaxy clusters, 
time-dependent terms A ~ ~ O{co'^) are important, 
where lo = 0{\k\a) is the inverse of the timescale to cross 
a system of size |k|^^ by stars of velocity dispersion a in 
unit of the speed of light. There can also be diffusion on 
small scale due to a pressure-like term V^A = — |kpA. 

The scalar field equation of motion becomes 



[9i 



(1 



cl)r]dt 



where 77 ^ is a damping time scale due to coupling of 
with the M field (cf. Appendix), and the diffusion 



term — c^V^ 



cf„|k|^ can be neglected if = 0. The 



scalar field ip then follows the equation of a damped har- 



monic oscillator with a damping rate (1 



)?7 and a 



slightly non-linear restoring force ^iPip and an exter- 
nal force ^ v'^^pM ^ |V<i>p(27V^(^M)- Assuming that the 
correction due to Hubble expansion ^ is negligible for a 
very small 6', the scalar field ip eventually approaches the 
MOND-likc static solution ipM: thanks to the damping 
term with a timescale ?7~^, which kills any history de- 
pendence. Rapid oscillations will likely keep the fluid's 
time-averaged property close to MOND-like solution as 
well. 

We estimate the oscillation time scale 



(27V)-i(AoF"A)-V2|, , 



10' 

^ • 300yr, 
(41) 

which is about 10' years if ^ 10. Here we assume 
Pm = VA = 0(1) = F" for systems of mild grav- 
ity 10^®cm/s^, e.g., clusters; for systems of stronger 
gravity, the time scale is perhaps longer). In the pro- 
cess of damping there will be a correction to MOND fiM 
function by the q term, heuristically, 1 — = 1 — ip\j if 



to respond to the formation of galaxies. Rather Lp would 
lack behind, might remain close to its cosmological aver- 
age: 

~ (43) 

which would mean a boost of the gravity of the baryon 

by a constant factor (1 — 'p^)~^ everywhere. 

7. GENERIC PROPERTIES OF DARK FLUID 

It is still uncertain whether the time-dependent cor- 
rection and a possible diffusion term are enough to help 
MOND to explain the Bullet Clusters (Angus et al. 2007, 
Angus & McGaugh 2008). However, it seems robust that 
the Dark Fluid - described by the field = I&°-p{\) 
is generally out of phase from the baryonic fluid. There 
are two types of deviations from MOND in general: 

• The Dark Fluid has a natural oscillation on time 
scales of v~^, which can be damped on a cross- 
ing timescale unless the external forcing is in res- 
onance. A very fast damping would mean an al- 
most instantaneous relation between gravity and 
the scalar field 1 — as the jiM in classical 

(39) Bekenstein-Milgrom (1984) modified gravity inter- 
X pretation of MOND. A slow damping would mean 
' a history dependent relation, reminiscent of Mil- 
grom's modified inertia interpretation of MOND: 
the dark fluid adds a dynamically-varying ineria 
around the baryons which it surrounds. A possi- 
ble test could be in galaxies with rotating bar(s), 
where there could be a phase lag between the bar 
and the effective Dark Matter (Debattista & Sell- 
wood 1998). This has intriguing consequences to 
the bar's pattern speed because of non-trivial cor- 
rections to the MOND pictures of dynamical fric- 
tion (Ciotti & Binney 2004, Nipoti et al. 2008, 
Tiret & Combes 2008); the properties of the dark 
fluid is in bewteen that of real particle dark halo 
and that naively expected from MOND. 

• The Dark Fluid has a pressure, controled by a prop- 
agation speed c^, where the speed of light is unity 
here, and the Dark Fluid can be made Cold by 



0, or Hot by c; 



1, or Superluminal by 



> 1. The ip would no longer be a function of 
the local gravity at r (as in MOND), rather it is a 
weighted average of a volume of all points ri by 
a Yukawa- type screening function exp(— '^^'"J"'"^ ), 
where 

108 



Ao 2N^Ao' 



(42) 



where q = [df — c^V^ + (1 — c'^)ridt\ is an operator. In 
tidally acting systems the value for ip will oscillate be- 
tween its pre-merging value and its equilibrium value. 

Models with a small N would not give MOND. E.g., if 
iV = 1 - 10, TTZ^-i ~ (100 - 10) Gyrs, then the universe 
would be too young dynamically to have a precise MOND 
effect in galaxies because p would not have enough time 

^ Considering the expansion of the universe would introduce a 
correction term Jb' in the force, where J = (3H^ + 2if?j) . 



Screening Length = Cy^u c,^-^ x 3001ight years. 

(44) 

Note this spatial correction to MOND can exist 
even in static systems; even a small pressure term 
with ^ might smooth out MOND effects 
on small scale structures (wide binaries, star clus- 
ters, dwarf galaxies), where the wavenumber |kp is 
much bigger than in galaxy clusters. The screening 
length can be set at ^ 100 pc for either a iV ^ 10^, 
~ 3 X 10^ km/s Hot Dark Fluid or a iV ~ 10'' 
and 30 km/s Cold Dark Fluid. This scale 

lOOpc is a scale dividing dense star clusters and 
fluffy dwarf galaxies. Observationally dark mat- 
ter effects are only seen in the universe on scales 



larger than lOOpc. It has been chaUenging for 
MOND to explain this observed scale (Zhao 2005, 
Sanchez- Salcedo & Hernandez 2007, Baumgardt et 
al. 2005). 

In conclusion, we find a framework of Dark Fluid the- 
ories where MOND corresponds a special choice of po- 
tentials or mass for the vector field. The Dark Fluid 
can run Cold or Hot depending on the sound speed 
(which could even be a running function of the vector 
field). These theories degenerate into scalar field theo- 
ries for Dark Energy effects in the Hubble expansion. It 
is possible to create an exact w = —1 Dark Energy effect 

(and a Dark Matter effect for b^^ = > 1). The 



scale uq — \/Aq in MOND in equilibrium spiral galaxies 
derives its physics from the amplitude of the dark energy 
Aq. MOND or Dark Matter effects are hence indications 
of a non-uniform Dark Energy fluid described generally 
by a vector field Z°'. For non-equilibrium systems like 
the Bullet Clusters or galaxies with satellites, the prop- 
erties of the Dark Fluid do not follow exactly the usual 
expectations of MOND or Cold/Hot Dark Matter, but 
(not so surprisingly) in between. 
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the figures. 



APPENDIX 

APPENDIX: THE VECTOR FIELD EQUATION AND THE DAMPING RATE 77 

Consider the vector field perturbation Xj — Yj ^ 1 in the compressional spin-0 mode with a potential Y. Define 
Co = — = ^TT-, we note the vector field EOM in the raw form is 



2A(y, + A,)\ - d, [2C2X {Y,, + e)] = 



(Al) 



If If 

for the index i = x,y,z, where n — 1~ c^, and Ai = u^VaWi and 9 = Vaw", and m° is the unit four-velocity vector 
locally. Apply the approximation A is a very small constant (strong gravity), so that Ai 



-<i> , and d,9 



-S^" ,:, and 



neglect the term — because — |aei| ^ 1 for perturbations, and |(^| <C inside the causual horizon. Replace 



if : 



\/\ and apply di to both sides, we get 



didt 



2A(r, - $,,) 



did. 



2C2X{Y,j, - 3*) 



4iV2A 



A. 



O 



IkPwA. 



Furthermore, we neglect the spatial and temporal variations of A and C2, factor out 2A, and define 



E 



a;" 



3 _ 



then the approximate equation for rj is obtained: 



[dtdt - CaV^] 77 ^ 



(A2) 



(A3) 



(A4) 



where n 



(? and C2 plays the role of sound speed squared. Replace ^ 



-<I>, and replace the partial derivative 



di with the wave vector k, and dt with the orbital frequency we get 

(l + 3C2)|k|2$ I + 3C2 



;2 + C2|k|2 l-C2a-2 



X 



0( [orbit crossing time] 



(A5) 



This estimation is good for fairly hot system a ^ 1/N ^ 30km/s if A^ 10*, and a ^ -^Z | C2 1 . This is an over-estimation 
if the vector field has a large (relativistic) sound speed due to a not-so-small •\/|C2|; this is an under-estimation if 
the vector field has a small sound speed y^jC^ in resonance with the stellar velocity a or if the system is very cold 
(7 < 1/A'^. There will also be corrections of order H/u in an expanding universe. 
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